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I. 



INTRODUCTION 



In any inventory model, one of the main objectives is 
to avoid both overstock and out-of-stock situations. The 
latter is especially critical in military supply. In lot 
size - reorder point models with stochastic demands, a 
decision has to be made as to when to reorder and how much 
the reorder quantity should be. The goal, then, is to 
achieve in some sense an optimal reorder level r and an 
optimal reorder quantity Q. One way to accomplish this goal 
is to define a cost function K which then acts as an objective 
function to be minimized by choice of Q and r as variables. 

Hadley and Whitin [Ref. 1] present several lot size - 
reorder point models. The one of interest in this paper is 
called a backorder model wherein, if demands occur when the 
system is out of stock, they will be backordered with some 
penalty cost, to be called the backorder cost. Customers 
then must wait until all orders are delivered to the inventory 
manager. Associated with delivery is a delay in time called 
lead time. In this paper we consider only the case where 
lead-time demand is assumed to be normally distributed. 

Mathematically, the backorder model of interest here has 
the objective function, 

K (Q,r) = + IC [^ + r- y] + ^-[ (y-r) $ + a< j> 1 (D 
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where 



A is the set-up cost (cost of placing an order) 
it is the backorder cost per unit backordered 
X is the average annual demand 
I is the inventory carrying charge 
C is the unit cost of the item 

y and a are the mean and standard deviation of the 

normally distributed lead-time demand random 
variable 

<J>(z) is the standard normal density function 

<Mz) is its complementary cumulative distribution function, 

/ 

oo 

i.e. , $ (z ) - f <j> (x) dx . 

z 

K(Q,r) is called the average annual cost and is composed 
of three different types of cost: the ordering cost, the 

cost of carrying inventory and the cost of backorders. These 
costs are represented, respectively, by the first, second and 
third term in the expression for K(Q,r) in (1). 

For ease in the algebra, (1) will be rewritten as 

K(Q,z) = + IC[| + az] + ^ [<Mz) - z<Mz) ] (2) 

where z = . 

Hadley and Whitin [Ref. 1, pp. 165 ff.] point out that 
the model just defined is appropriate when the expected number 
of backorders is negligible. They then propose a solution to 
the optimization problem by claiming that K is convex in Q 



7 



and r (or z) and therefore any solution obtained by setting 
the partial derivatives of K equal to zero will determine a 
global minimum. They provide an iterative scheme for such 
a solution along with the assertion that the solution, when 
it exists, is unique. Finally, they discuss the fact that 
a solution will not always exist, pointing out that this 
"anomaly" arises because a backorders term was omitted in the 
expression for the carrying cost. Since the model is to apply 
only when such a term is negligible, they assert that this 
situation will never occur when high backorder costs are 
involved without being explicit as to what "high" means. 

To verify these results, the reader is directed to solve 
a series of exercises, the first of which is to show that 
the function involving the backorder cost (the third term) 
in (1) , say J(Q,r) , is convex in Q and r. It would then 
follow immediately that K is convex in Q and r. But it was 
pointed out as early as 1964 by Veinott [Ref. 2] that J is 
not convex. In 1969, Brooks and Lu [Ref. 3] addressed the 
same problem and derived a general result which, when applied 
to the normal case, shows that J is convex for all r _> p , 
the mean lead-time demand. 

However, these results still leave the main issue 
unresolved. In the first place, solutions can be obtained 
where r < p.'*' More importantly, the real issue is whether 

■'"The difference r-p is sometimes referred to as a 
s afety level . In that case, a negative safety level simply 
means it is more economical to tolerate a backorder position 
on the average than to protect against that position. 
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or not K is convex, and, in any case, how its minimum value 
is characterized. The purpose of this thesis is to clarify 
this issue and characterize the solutions to the problem in 
a more explicit manner. 

Section II will show that in fact K is not convex and, 

in addition, will exhibit an example in which the first 

d K dK 

partial derivatives and vanish at two different 

points. Thus, the solution is not unique as claimed by 
Hadley and Whitin. Section III will deal with the character- 
ization of the solutions and will show that under suitable 
conditions there are always two distinct solutions, one of 
them being a minim-urn and the other a saddle point. Conditions 
for the existence of a solution will be given in terms of the 
set-up cost A. Section IV will discuss the application of 
the results in Section III to practical inventory problems 
along with the iterative scheme given by Hadley and Whitin to 
solve for a solution. 
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II. COUNTEREXAMPLES 



As pointed out in the introduction, the backorders term 
in the cost function is not convex as asserted by Hadley 
and Whitin [Ref. 1, Problem 4-6, p. 221]. Even so, it might 
still be the case that K is convex as asserted by the authors. 
That this is not so will be demonstrated in part A. On the 
other hand, even though K is not convex, it might still be 
the case that the partial derivatives vanish at a unique 
point and that point determines a minimum. That this is not 
so will be shown in part B and a characterization of solutions 
will be given in the next section. 



A. CONVEXITY CONDITION FOR K 



The function K is convex if and only if its Hessian 
matrix ft is positive semi-definite [Ref. 4]. Since K is a 
function of two variables, the matrix ft will be positive 
semi-definite when its diagonal terms and its determinant 
are nonnegative. 

For the case at hand, using (2) , the Hessian matrix for K 
is given by 



ft = 



^ ( A + Trn (z) ) 
Q J 



itX<I> (z) 

Q 2 



TT A 4> ( z ) 



Q“ 



it A 4> ( z ) 
cQ 
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where n(z) = 0<J)(z) - oz<Hz) 

or tT ( z ) = au(z) (3) 

where u(z)=<J>(z)-z$(z) . (4) 

In Appendix A it is shown that u(z) > 0 for all z . 
Hence, for nonnegative set-up costs A, it is clear that the 
diagonal terms of Q are positive. Thus K is convex if and 
only if det ft j> 0 . 

Now, 

det n = (A + „n(z>> - 

oQ Q 4 

and thus det ft _> 0 if and only if 

„ a 2 x/ X 2,2.2, , 

? JL ( A + tn( 2)) - * x ♦ (z > > o 
oQ Q 

4 . . 

Multiplying the above inequality by Q > 0 and dividing 
2 

by tt X > 0 gives 

1 L $ . J Z ) (a + tt rf ( z ) ) - tt$2(z) o 

O A 9 9 

or — (Hz) + 2tt<H (z) - 2ttz(J) (z) $ (z) - tt<)> (z) 0 (5) 

Defining V(z) to be equal to the left-hand side of (5) , then 
the following summarizes convexity. 
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RESULT: 



The function K is convex if and only if det ft > 0 

or equivalently if and only if V(z) >_ 0 . 

REMARK : This result holds for any Q > 0 and it is shown 
in Appendix A that V(z) > 0 whenever z >_ 0 . 

A counterexample to the convexity of K is then attainable 
from an example provided by Hadley and Whitin [Ref. 1, p. 173] 

There the parameters chosen are tt = 2000, X = 1600 , C = 50 

I = .20 , A = 4000 , y = 750 , o = 50 . For such a case 
V(z) becomes 

V(z) = 1 6 0 (J> ( z ) + 4000(J) 2 (z) - 4000z<J> (z) 4> (z) - 20004> 2 (z) 

and for z = -1 (equivalently, r = 700) 

V(-l) = 38.72 + 234.256 + 814.3784 - 1415.57138 
= -328.21698 < 0 



Thus, K cannot be convex. 

B. COUNTEREXAMPLE TO UNIQUENESS OF THE SOLUTION 
DEFINITION. (Q,z) will be called a solution whenever Q > 0 
and z satisfy the following equations: 



0 



3K 

3Q 




ttXti (z) 





nXau ( z) 




( 6 ) 



0 



9K 

3z 



(IC - 



7TX<i> (z) 

Q 



) O 



(7) 
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It is convenient to write (6) and (7) as: 



-if 



IC 



-f 



2A(A + iTri(z)) _ ,I2A(A + 7Tau(z)) 



IC 



and 



4>(z) = 



QIC 



— 



( 8 ) 



(9) 



COUNTEREXAMPLE . Let tt = A = o = 1 , I = .1 , C = 2.42 

12 

and y = 3. Let A = ~ 2 ' ( ' i ' ) ~ ( ' 2 '~ 4 2 ' )" ^ “ u(.05) = .1018 , 

where u(z) is defined in (4) . Then there are two distinct 
solutions 

< ^(z 1 ) 

1) (Q 1 ,z 1 ) = ( J7XT(2T42) r z l ) w here z ][ = .05 

$(z 2 ) 

2) (Q 2 ,z 2 ) = ( y i) (2 ~ 42 ) ' z 2 ^ where -1.81 < z 2 < -1.80 

To see that (Q^' z ^) is a solution. 




Z 1 



•j r $ (.05) , AE -^ i 

2 ( ,1) (2.42) UUU ^ J + (.1) (2.42) _ (1) (1) (l)u(.05) 

4> 2 ( .05) 2 4> 2 ( .05) 

( .1) 2 (2.42) 2 ( .1) 2 (2.42) 2 

( .1) (2.42) A ( .1) 2 (2.42) 2 u( .05) A (.1) (2.42) (.l) 2 (2.42) 2 u( .05) 

2 f o f 2 

4> ( .05) Z $ (.05) 



= 0 . Thus, equation (6) is satisfied. 
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a k 
TZ 



= (.1) (2.42) - ~ ) ff^y - ° - 5 ' = (.1) (2.42) - (.1) (2.42) = 0 • 
Q 1 ,Z 1 (.1) (2.42) 

Hence equation (7) is satisfied. Therefore, (Q^,z^) is a 
solution . 

1 2 

Now, let S (z ) = 2 ( l) (2 ~ 42) ^ ^ • Since S is 

clearly a continuous function and 

S (-1 . 81) £ .0997 < A < .1060 <_ S(-1.80) 

there must exist some z 2 such that -1.81 < z 2 < -1.80 
and S(z 2 ) = A. The verification for (Q 2 ,z 2 ) is done in a 
similar manner. 



a K 



= _ (-D (2.42) (.1) 2 (2.42) 2 u( z 2) (.1) (2.42) _ ( .1) 2 (2 .42) 2 u( z 2) 

2 4> 2 ( z 2 ) 2 4> 2 ( z 2 ) 

=0 . Thus, equation (6) is satisfied. 



2 

* (Zo) 

1 [ 2(.j) (2.42) ~ u(z 2 )] 
$ 2 ( z 2) 



(.1) (2.42) (1) (1) (l)u( z 2) 



4> 2 ( z 2) 



( .1) 2 (2.42) 2 



( .1) 2 (2 .42) 2 



3K 

8z 



Q 2 ,Z 2 



(.1) (2.42) 



= - ( *D (2.42) =0 

(.1) (2.42) 



and hence equation (7) is satisfied. Therefore, (Q 2 ,z 2 ) is 
also a solution. 
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To illustrate these results, a graph similar to the one 
given on p. 171 of Ref. 1 is reproduced in Fig. 2-1 for this 
special case. In this graph r = az + jj is expressed as a 
function of Q for each of Equations (8) and (9) . It shows 
that the two curves intersect at more than one point. 
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1 2 

The value A = ' 2 ' ( ' l) ' (2 42) ^ (- 05 ) “ u(.05) was by no 
means chosen by accident. As it will be seen in the next 
section, the set-up cost A plays an important role in the 
characterization of the solutions, and it will be shown that, 
under appropriate conditions of the other parameters, there 
are always two distinct solutions, one of which is a minimum. 
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III. CHARACTERIZATION OF MINIMUM SOLUTION 

Let a set of parameters be given. To simplify matters, 
define 



A series of lemmas are convenient to derive the main result 
in a theorem. 

LEMMA 1. Let (Q*,z*) be a solution. Then, 

(a) A = S (z*) 

(b) V(z*) > 0 if and only if <(>(z*) > and V(z*) = 0 

if and only if <{>(z*) = . 

T TT A 

Proof (a). By definition of a solution, Q* and z* must 
satisfy (8) and (9) of Section III. Thus, 




( 10 ) 



2 it (A+ttctu (z*) ) 
IC 



(Q*) 



2 _ (ttA) 2 $ 2 (z*) 
(IC) 2 



or 



2AA + 27TAau(z*) = $ 2 (z*) 



2AA = $ 2 (z*) - 2 ttActu (z*) 



A = 4> 2 (z*) - TTCTU ( z* ) = S ( z*) 
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Proof (b) . 



s> a ? 

V(z*) = — <p (Z*) + 2TT<p(z*) ( <f> ( z*) -z*$ ( z*)) - TT$ (z*) 



2A 2 

= =£■ <p(z*) + 2tt4>(z*)u(z*) - TT<r(z*) 



= 2 <f> ( Z * ) (| + TTu ( Z * ) ) - TT$ 2 (z*) 



r 7T 2 X ,2 



= 24 >(z*) (' 2x"cq ^ ( z *) -Tru ( z *) +Tru (z*) ) - tt$ (z*) from part (a) 



^ 4»(z*) $> 2 (z*) - tt $> 2 (z*) 



^ $ 2 ( z *) t4>(z*) - 



Since 0^(z*) > 0, the result stated in (b) 

immediately . 



f OllOV7S 



Q.E.D. 



It is not always true that there exists a solution 
for any given set of parameters. The following lemma discusses 
conditions for which K cannot be minimized. 



LEMMA 2. Given A >_ 0 and tt, X, I, C, o all positive, a 
necessary condition for the function K to be minimized is 
that < <fr(0) . 

Proof. The proof is divided into two cases. 

ICo 

Case (i) : Suppose — > <J>(0) and K is minimized at (Q*,z*). 

~ TT A 

Then the partial derivatives must vanish at (Q*,z*) . But 
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V ( z*) 



= 15? < 0 

since $(z*) jf. <f> C 0 ) and so (Q*,z*) is a saddle point. This 
is a contradiction. 

Case (ii) ; Suppose yy = 4> ( 0 ) and K is minimized at (Q*,z*) 
Again, the partial derivatives must vanish at (Q*,z*). Now 
if z* / 0 then V(z*) < 0 and hence (Q*,z*) is a saddle point 
Thus it must be the case that z* = 0. But then, by Lemma 1, 
part (a) , 



TT^ X 2 

A = ^ <*> ( 0 ) - TTCU(O) 



ETE' - *04,(0) 



TTO 

8 



irX 

ICa 



- ira$(0) 



no 1 

8 ’ 4 .( 0 ) 



7TO$ (0) 



= 7TC ( 



1 - 84> 2 ( 0 ) 

$70] 



) 



2 

Since 1-8$ (0) < 0, A < 0 contradicting the fact that 
A j> 0. 

ICa 

Thus K cannot be minimized unless yy < $(0) . 

Q.E.D. 
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ira 

LEMMA 3. Suppose — y- < 4>(0) . Let z^ > 0 and -z. < 0 

TT A O O 

. ICQ 

be the unique points such that <J>(-z ) = <j>(z ) = — r- . 

Then, with S(z) defined in (10), 

(a) S has a minimum at z Q where A^ = S(z Q ) < 0 

(b) S has a maximum at -z Q where A Q = S(-z Q ) _> 0 

(c) S is strictly concave increasing on the interval (-°°,-z ) 

ir^ X 2 

Proof. S ( z) = $ tz) - iTau(z) 



S ' ( z) = - 



TT 2 A 
IC 



<J>(z) <Mz) - TTCTu' (z) 



TT 2 X 
IC 



<J>(z)$(z) - rcr(-$(z)) (see Appendix A) 



= ttct$(z) - 



TT 2 X 
“IC 



<t> (z) $ (z) 



9 

tt X . . . r ICa , . . . 

= Tc" t<z) [ inr - * (z >i 

. 2 , 

Since $(z) > 0 , the first derivative S'(z) vanishes 

only at -z and z Q . Also, S' (z) < 0 for -z Q < z < z Q and 
S' (z) > 0 for z < -z or z > z Q ; that is, S is decreasing 
in the interval (-z 0 ,z Q ) and increasing in the intervals 
(_co,_z o ) and ( z q , °°) . 
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2 2 

s " (z) = - t(z) - 4>(z)) + <D(z) (z<J>(z) ) 

= -ira<j> (z) + ? 2 - A ^ (z) + ll.A_z_^z)^(z) 

2 

= — ( Cf> ( Z ) + z<±> ( z) ) - Tra<J>(z) 

= ira<J>(z) [~ (<p (z) + z*(z) ) - 1] 

Thus, at z Q , 



S"(z 0 ) = 






<^f + 



z o* (z o }) - 



1] 



= ira <f> (z D ) [1 + 



ttA 

ICa 



z o $(z o ) 



1 ] 



TcT z o tf>(z o )0(z o ) > 0 



( 11 ) 



Therefore, S has a minimum at z . 

o 

To show that A^ = S(z Q ) < 0 , suppose, to the contrary, 
that S ( z Q ) 0 . Then, since S'(z) > 0 for z > z Q , 

S(z) > S(z Q ) 0 for all z > z Q and hence lim S(z) > 0. 

Z->oo 



But lim S(z) = lim $ 2 (z) - vo lim u(z) = 0 (see 

Z-*°° Z->oo 1 C z^-oo 

Appendix A) . Thus, A^ = S(z Q ) < 0 . This together with 
(11) proves part (a) . 
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Now at -z 



o ' 

s"(- Zo > - ™*<-z 0 )[rar (^f - 2 o^- z o»> - n 

“ TO0(-z o ) [1 - z 0 4(-z 0 ) - 1] 

A 

= - -Jq- z 0 *(-z 0 )M- Zo ) < 0 . (12) 

Therefore, S has a maximum at -z . 

To show that A q = S(-z Q ) 0 , again suppose that 
S(-z ) < 0. Since lim S(z) = 0 , there exists some z' 

2-^-co 

such that S(z') > S (-z Q ) , for otherwise S(z) <_ S (-z ) for 

all z implies that lim S(z) _< S (-z Q ) < 0 . But the existence 

z -)-co 

of such z' contradicts the fact that S is maximized at -z . 
Thus, A = S (— z ) > 0 . This result coupled with (12) proves 

part (b) . 

Part (c) can now be easily proved by noting that for 
z < -z Q < 0 , 

<t>(z) + z$(z) < <t> ( z) < $(“ z 0 ) = 
hence -j—- ( <J> ( z ) + z$(z)) < 1 
or (<j>(z) + z$(z)) - 1 < 0 

and therefore S" (z) = ttct (z) (<|>(z) + z$(z)) - 1] < 0 

Thus, S is strictly concave increasing on the interval (-»,-z ) 
The proof of Lemma 3 is now complete. 

Q.E.D. 
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REMARK : The minimum and maximum in Lemma 2 are clearly unique. 



Moreover, lim S(z) = (see Appendix A). A graph of S(z) 

z-*~ 00 

corresponding to the example of Section II is given in Fig. 3-1. 




The following lemma shows the importance of A in the 
characterization of a solution. 

ICa 

LEMMA 4. Suppose — r~ < cf>(0) and A > 0 is given. Then 

TT A — 

(Q*,z*) is a solution if and only if A = S(z*). 

Proof. The necessary part has been proved by Lemma 1, part 
(a) . The sufficient condition will now be proved. Suppose 

A = S ( z*) = ^ <j> 2 ( 2 *) - ttou ( z*) . Let Q* = . 

Then, it can be easily shown by substitution of A and Q* 
into equations (6) and (7) that (Q*,z*) is a solution. 

Q.E.D. 
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REMARK : It should be observed immediately from Lemmas 4 and 3 



that a necessary condition for a solution to exist is that 
0 <_ A _< A q . The following theorem summarizes all the results. 



T P rr 

THEOREM. Given tt, X, I, C, a all positive and — < <f>(0) 

7T A 

and A q > 0 , then for 0 < A < A q there exist exactly two 
solutions (Q^*,z^*) and (Q 2 * r z 2 *^ such that z^* > Z 2 * ; 
(Q^*,z^*) is a minimum and ^2* ,z 2*^ a sa< 3dle P°int for K . 



Proof. Since S(z) is strictly concave increasing in the 



interval (-<»,-z ) , and strictly decreasing in the interval 

(-z 0 ,z Q ) with lim S(z) = , S (-z Q ) = A q > 0 and 

z -*— 00 

S(z Q ) = A^ < 0 , the function S(z) must have two roots z^ 

and z 0 where z-, < -z_ and -z_ < z- < z^ . Thus, for any 

2 1 o o 2 o 



A such that 0 £ A < A Q , 



there are two distinct values 



z^* and Z 2 * such that z^ < z^* < -z Q < z-^* < Z 2 and 



A = S(z, *) = S(z~*). By merely defining Q 1 



* = 



7rX$ (z^* ) 



IC 



and Q 2 * = 



ttX 4> (z 2 *) 



IC , it follows by Lemma 4 that (Q^*,z^*) 

z_ implies 



and (Q2* /Z 2*^ are s °l ut i° ns * But, -z Q < z^* < _ c 

that 4> (z 1 *) > arid hence by Lemma 1, part (b) , V(z^*) > 0. 

Thus, (Q^*,z^*) minimizes K. Similarly, Z 2 * < -z Q implies 

ICa 

that 6( z~*) < — and hence by the same lemma V(z 0 *) < 0 . 

2 7T X 2 

Thus, (Q~* ,z 0 *) is a saddle point. 

* Q.E.D. 

REMARK: For A = A Q = S(-z Q ) , the only solution possible is 



determined by z 



* = - . » 



But then, 4>(z*) = 



ICc 

TT X 



and thus 



V(z*) =0, so that no conclusion can be drawn. 
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IV. COMPUTATIONAL ASPECTS 



Before attempting to solve a given problem involving a 
set of parameters, some simple tests should be performed to 
see whether or not there is any solution . 

1) The first test is ^2- < <f> ( 0 ) . 

2) If this test is passed, then find -z and A = S(-z ) . 

o o o 

Two cases could then occur. 

Case (i) : A q = 0 . The only solution then is -z and V(-z ) =0 

and it has been pointed out that no conclusion may then be 
drawn . 

Case (ii) : A q > 0 . If 0 < A < A q , then there exist exactly 

two solutions (Q^*,z^*) and ($2* w ith z 2* < " z 0 < z l* * 

As shown, (Q^*/Z^*) is a minimum and (C>2* ,z 2*^ a sa £dle 

point. If A = A , then again the only solution is -z and 
no conclusion can be drawn as in case (i) . 

Hadley and Whitin give an iterative scheme for solving 
for (Q*,r*) in Ref. 1. This scheme seems to converge always 
to the minimum value and never to the saddle point. But there 
is a simple explanation for this. The scheme begins by choosing 

initially a Wilson Q w = for the value of Q to be 

used in Eq. (9) (see Fig. 2-1) and Q w < Q^*. The iteration 
then proceeds to choose a value of r corresponding to from 
Eq. (8) . With that value of r a new value of Q is selected 
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from Eq. ( 9 ) and it will always be the case that Q < Q^* 
in that selection. Proceeding in this manner, the pairs 
(Q,r) converge to the solution (Q^*,r-^*) and cannot yield 
(Q2*,^2*). Thus, when a minimum exists, this computational 
scheme may be relied upon. 

However, in any case of doubt, a simple test is furnished 

ICo 

by checking whether or not > yy for the solution. 

It has been pointed out that the boundary case A = A Q 
causes problems in that no conclusion may then be drawn. This 
presents no difficulty for practical uses, however. After 
all, an optimal solution may be given for any A arbitrarily 
close to A q . In an application it would be a rare circumstance 
that produced a set of parameters in which A (typically 
measured in dollars and cents) would be exactly A Q . 
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V. CONCLUSION 



While the results of this thesis only apply to the case 
of normal lead-time demand, they are of more than mere 
academic interest. In fact, that is a common assumption made 
by Hadley and Whitin along with users of the formulas and 
schemes given by them. And it is a happy coincidence that 
the algorithms given are indeed valid even though the 
assumptions under which they were derived were replete with 
errors . 

Nevertheless, it is a recommendations of this thesis that 
further study into other types of lead-time distributions be 
made. One of the key features of the results derived here 
was the examination of the Hessian matrix with the observation 
that its positivity did not depend on Q. This will still be 
true with a density other than normal and so much of the same 
general technique ought to apply. 

One of the features of this thesis has been to amplify 
and clarify completely the minimum solution. Conditions are 
given for existence, and formulas are given along with, simple test 
criteria when existence is assured. A by-product of the 
results is the fact that negative safety levels do indeed 
occur as the unique optimal solution. In fact, a simple 
examination of a graph like S(z) will show immediately and 
simplify those values of the set-up cost A for which such 
solutions are valid. Since the solution in such cases is 
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unique, certainly negative safety levels should not be ignored 
or tampered with (truncation to y for example) in an application. 
It should also be noted that the solutions are characterized 
independently of y. 

As a final recommendation, this thesis was restricted to 
the so-called backorders model where a customer must wait for 
delivery when the inventory position is out of stock. There 
are, however, other models such as the "lost sale case" 
discussed by Hadley and Whitin to apply in this situation. 

An analysis similar to that applied here is surely called for 
so that resultant solutions are validated and characterized. 

For the authors continually appeal to convexity for drawing 
conclusions in those cases. 
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APPENDIX A. VARIATION OF V(z) 



As in Section II, 



ft 



—■ (A + 1T1-) (z) ) 
Q 



pA$ ( z) 



ttA$ (z) 

.2 



Q' 



Q 



ttA < f> (z) 

OQ 



The diagonal term — is certainly positive. 
W(Z) = o(<J>(z) - z$(z)) = ctu(z) as defined in (4) 

u' (z) = — z 4> ( z ) - (4>(z) - zcf>(z)) 

= -<Mz) < 0 for all z 

u" (z) = <t>(z) > 0 for all z 



Clearly, lim u(z) = lim ( 4> ( z ) - z$(z)) = 0 + 

2-v- CO Z->-oo 



and 



lim u(z) = lim (<|>(z) - z$(z)) = - lim 

Z-voo £-voo £->cn 



= - lim 

Z -M» 



Z 



1 

$ (z) 



= -J-00 



z$ (z) 
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Applying L'Hospital's Rule repeatedly, 



lim u (z) 
z-*-°° 



lim 

z~*~ m 



1 

( Z ) 

Vui 



lim 

Z-+°° 



$ 2 (z) 
<t> (z) 



lim 

Z^-oo 



-2$(z) <fr(z) 
-Z<$>(z) 



lim 

z ->oo 



2 $ (z) 
z 



0 . 



Thus, u(z) is a strictly convex decreasing function 
approaching zero as z + m and hence u(z) > 0 for all z. 
A sketch of u(z) is given in Fig. A-l. 




Fig. A-l 



Since u(z) > 0 for all z, n(z) 
and therefore the diagonal term 
all z since A > 0 . 



2 X 



au(z) > 0 for 
(A + irn(z)) > 0 



all z 
for 
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2A 2 0 

V(z) = — <J>(z) + 2 tt (z) - 2ttz<|> (z) $ (z) - 7T<J> (z) 



Define T(z) = 2 <f* 2 ( z ) - 2z$ (z) 4> (z) - $ 2 (z) , then 



V (z) = 4 (z) + 7rT (z) 



T'(z) = -4z<f> 2 (z) + 2z 2 <{> (z) 4> (z) + 2z<f> 2 (z) 



= 2z 2 (j) (z) <J> ( z) - 2z<J> 2 (z) 



2 z 4> (z) (z<J> (z) - <}> (z) ) 
-2z4 (z) u ( z) 



T* 1 


(z) = 


= — 2 4> ( 2 


Ou(z) + 


2z 2 (J) (z) 


Thus , 










V' 


(z) = 


2A 

0 


z (z ) + 


ttT ' ( z ) 






£51° 

I 

II 


z<J>(z) - 


2ttz<J) ( z) 






= -2z<J>(z) ( ^ + 


iru (z) ) 


Since 


u (z) 


' > 0 


for all 


z 




j 


r ' ° 


for 


z < 0 


V' 


(Z) / 


' = 0 


for 


z = 0 






1 < o 


for 


z > 0 



from (A. 1) 



(A.l) 
(A. 2) 
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V" ( z) » 



2A ? A 9 

- — <f>(z) + — Z ( z ) + ttT" (z) 



At z = 0 , 



2A 

V" (0) = - ~ (0) + ttT" (0) 



2A 



<£(0) ~ 2tt(J) (0) u (0) from (A. 2) 



Thus V" (0) < 0 and hence V(z) is maximized at z = 0 and 
V(0) = ~ <M0) + ttT(O) = ^ (.3989) + tt(.0683) > 0 . 



lim V(z) = 0 + u lim T(z) = 0 

2^+°° 2^+oo 

lim V (z) = 0 + t\ lim T(z) = tt (-1 ) = - it . 

2->— oo — oo 

A sketch of V(z) is given in Fig. A-2. 




Since V(z) attains a maximum with positive value at z = 0 and 

is decreasing for z > 0 and lim V(z) = 0 , V(z) > 0 for all 

z-h-°° 

z > 0 . 
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